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GOTTINGEN 



Zusammenfassung 


ImRahmen des Forsch\ingsprojektes wurde 'ein System fur ■die Aus'wer- ' 
tung von Interferenzmustern entwickelt. Interferogramme von klassisch.en 
und holografischen Interferometern werden mit Hilfe eines Fernsehabtast- 
verfahrens digitalisiert und in einenRechner eingegeben. Nach einer von 
der Qualitat der Vorlagen abhangigen Bildvorverarbeitung ■wie zeitlicher 
und rS.umlicher Mittelung, Subtraktion der Untergrundintensitat oder z-wei- 
dimensionaler Filterung zur Elimination von Stormustern und einer Grau- 
■werttransformation werden die Flanken der Interferenzstreifen als Poly^ 
gonziige aus dem Interferenzbild extrahiert. Diese Darstellung wird nach 
der automatischen Numerierung der Linien und einer Beseitigung von 
eventuell durch Kontrast- oder Auflosungsmangel aufgetretenen Linien- 
unterbre,chungen in. eine Polynomdarstellung gebracht, die die Objektfunktion 
durch einen Satz von Koeffizienten an einem regelmafiigen Gitter beschreibt. 

Es werden Beispiele zur. Mittelung .und-Differenzbildung von Interferogrammen. 
so\yie ^zur Analyse von.Strukturen (Herste llungsfehlern), ,in . Xesthologrammen. 7 ^ 
von Kraftfahrzeugreifen angegeben. 



Summary 


During the research period of this project a system for the evaluation 
of interference patterns has been developed. For digitizing and proces- 
sing of the interferograms from classical and holografic interferometers 
a picture analysis system based upon a computer with a television digi- 
tizer has been installed. Depending on the quality of the interferograms, 
four different picture enhancement operations may be used: signal averaging, 
spatial smoothing, subtraction of the overlayed intensity function and the 
removal of distortion-patterns using a spatial filtering technique in the 
frequency spectrum of the interferograms. The extraction of fringe loci- 
from the digitized interferograms is performed by a floating-threshold 
method. The fringes are numbered using a special scheme after the re- 
moval of any fringe disconnexions which appeared if there was unsufficient 
contrast in the holograms. The reconstruction of the object function from 
the fringe field is achieved by a local polynomial least-squares approxi- 
mation with spline fit. Applications are given, demonstrating the averaging 
and subtraction of interferograms as well as the analysis of structures 
(defects of fabrication) in the holografic interferograms of car-tyres. 
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0, Einleitung 
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Die vorliegende Arbeit ist eine Zusammenfassung der dreijahrigen Tatig- 
keit an einem vom BMFT geforderten Projekt, das zum Ziel hatte, ein 
System zur automatischen Auswertvmg von Interferogrammen imd Holq-. 

; grammen zu entwickeln. 

Wahrend dieser Zeit wurden die einzelnen Verfahrensschritte wie Aufbau. 
des Rechnersystems und der Hardware, graph, Dateneingabe und -mani- 
pulation, Grauwertbildaufnahme und -digitalisierung, Interferenzstreifen- 
erkennung und Numerierung, sowie Erkennung und Beseitigung von even- 
i tuell aufgetretenen Auswertungsfehlern entwickelt und erprobt, Aufierdem 

wurde eine Methode der zweidimensionalen Polynomapproximation ent- 
wickelt um die ausgewerteten Interferogramme in eine numerisch geeignete 
\ Form zu transformieren, und sie der letztlich vorgesehenen Endauswertung I 

. zuzufuhren. . 

1 ? 
\ 

Die angesprochenen Verfahrensschritte sind ausfiihrlich bereits in den : 

Zwischenberichten I-V [1] zum Forschungsprojekt beschrieben worden, 
so da(3 sie hier nur noch relativ kurz besprochen werden sollen. 

Angewendet wurde das System fiir die Auswertung von Mach-Zehnder-Inter- 
ferogrammen von instationSren transsonischen Stromungen [6, 7], sowie 
fiir die Analyse von Doppelbelichtungshologrammen von Kraftfahrzeugreifen^[8] 
der Fa. Continental-Gummiwerke. 

f 

Fin betr achtlicher Teil der Zeit wurde fiir das Auffinden und die Analyse 
von Strukturen in den Reifentesthologrammen der Fa. Continental aufge- 
wendet, die Riickschliisse auf Reifenherstellungsfehler geben sollen. 
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Fur die Wahl des von uns verwendeten Fernsehabtastverfahrens zur Auf.- 
nahme der Interferogramme sprechen neben der Geschwindigkeit und des 
reichhaltigen Angebotes von Videokameras mit speziellen Aufnahmerohren 
auch die erzielbaren Auflosungen, die mit den zum Projektbeginn auf dem 
Markt befind lichen Bildspeichersystemen der in Fragekommenden Preis- 
klasse nicht erzielt werden konnten. 

Die mit guten Videokameras erzielbaren Auflosungen von 500 - 600 Linien 
pro Bildzeile sind fur die Auswertung der Interferogramme ausreichend - , 
wie auch eine Untersuchung an synthetischen Hologrammen gezeigt hat [2]. 
Die Grauwertauflosung ist im Vergleich mit Scanningphotometern bei , 
Videosystemen zwar nicht sehr grofi, reicht aber fiir die Analyse von 
morphologischenMerkmalen vdllig aus [3], .. 

Einige inzwischen auf dem Markt befind liche Bildspeichersysteme mit ver-. 
besserter Aufldsung konnten durchaus an dem System eingesetzt werden, 
zumal wenn sie durch geeignete Hardwareausriistung zu einer gewissen 
Bildvorverarbeitung (z.B. Bildmittelung) geeignet sind. Von uns wurde 
keine spezielle Hardware zur Bilddatenvorverarbeitung entwickelt, da 
der Schwerpunkt des Projektes auf der Auswertung der Interferogramme 
und der.en Weiterverarbeitung lag. 


1 . 


Beschreibung des GerS-tes 


Das System zur Auswertung von Interferogrammen besteht aus den Kom- ■ 
ponenten Rechnerhardware und -software, 'E in- lind Ausgabe graph. Daten, 
Grauwertbildeingabe, Grauwertbildausgabe, Filmprojektorsteuerung und 
Hologrammrekonstruktionsvorrichtung. Diese Komponenten werden im 
folgenden genauer beschrieben. Eine Konfigurationsskizze des Systems 
zeigt Abb. 1. 


1 . 0 Art der Vorlagen 

Die zur Zeit verwendeten Bildvorlagen sind mit Hochgeschwindigkeits- 
kameras auf 16 mm-Film aufgenommene Mach-Zehnder-Interferogramnie 
von instationSren transsonischen Strbmungen [6, 7] und auf Film regir 
strierte Doppelbelichtungshologramme von holografischen Kraftfahrzeug- 
reif enuntersuchungen der Fa. Continental Gummiwerke Hannover [8].. 

Die Abb. 3a und 4a zeigen typische Beispiele. Die Reifentesthologramme 
werden in einer Rekonstruktionsvorrichtung mit einem He-Ne-Laser re- 
konstruiert und mit einem Weitwinkelobjektiv auf die Aufnahmerohre der 
Videokamera abgebildet. Die Weitwinkeloptik ist erforderlich, weil 
jedes Hologramm einen ca. 1 00°-Ausschnitt der von innen aufgenommenen 
Reif enlaufflM.che zeigt. 

Grundsatzlich sind jedoch auch Vorlagen anderer Formate mit dem System 
auswertbar, da durch Wechseln des Kameraobjektivs der Abbildungsmafi- 
stab in einem grofien Bereich eingestellt werden kann. Ftir die Auswertung 
anderer Formate miifite lediglich die automatische Zuftihrung (Filmprojek- 
tor), die bei der Serienauswertung unumganglich ist, angepaBt werden. 
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Abb. 1: Konfigurationsskizze des Rechnersystems 























Die BildqualitS.t der Mach-Zehnder-Interferogramme zeichnet sich durch 
relativ gute Kontraste, sowie durch eine recht gleichformige Ausleuch- 
tung des Testobjektes aus, AuBerdem sind in den meisten Fallen keine 
StSrobjekte im Mefifeld. In einigen Fallen, in denen Markierungen im 
Mefifeld vorhanden sind, sind die jedenfalls ortsfest und konnen so leicht 
beriicksichtigt werden. 

Die Reifentesthologramme haben dagegen einen stark variierenden, zum 
Teil sehr geringen Kontrast, der zu den Bildrandern hin stark abfailt. 

Durch ungleichmaBige Ausleuchtung bei der Hologrammaufnahme und. 
Streuerscheinungen, die zu einem Interferenzmuster auf dem Hologramm- 
film fiihren, ergibt sich iiberdies eine starke Modulation der Untergrund- 
intensitat. Ferner zeichnen sich die Heizbalgabdriicke zum Teil sehr stark 
ab. Diese Heizbalgabdriicke sind schrag zur Reifenlaufrichtung verlaufen- 
de schmale Verdickungen, die bei der Reifenherstellung entstehen und 
nicht vermeidbar sind. Auf diesen Heizbalgabdriicken ist das Interferenz- 
muster z. T, bis zu einer halben Streifenbreite verschoben was u.U. zu 
Fehlinterpretationen bei der spateren Auswertung fiihren kann. Da durch 
das Aufnahmeprinzip das Objekt mit abgebildet wird, fiihren alle im MeB- 
bereich liegenden Objekte, wie die sog. Spreizer (die Stangen, die den 
Reifen auseinanderhalten) sowie ungleichmaBige Reflexion auf der Reifen- 
oberfiache zu zusatz lichen Stdrungen. Die Abb. 3b und 4b zeigen typische 
Kontrastver laufe der Interferogramme und Reifentesthologramme. 


1 . 1 Rechnerkonfiguration 

Die Rechnerkomponente besteht aus einem PDF 11/34 Rechner (Digital 
Equipment) der mit 25 6 Kilobyte (KB) MOS-Speicher ausgeriistet ist. Als 
externe Speichermedien stehen ein Floppy-Disk-Doppellaufwerk (je 240 KB) 
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sowie zwei Wechselplattenlaufwerke RLOl (je 5 MB) zur Verfiigung. 
Weiter ist der Rechner bestiickt mit einem 16 Bit Para lie linterface 
zumAnschluB eines Kurvendigitalisierers zur Eingabe graph. Daten, 
mit einem 16-Bit DIVIA- Interface zum Anschlufi des Bilddigitalisierers 
und mit einem seriellen Interface (V 24) zur Kommunikation mit einem 
Rechnerverbundsystem. Zur Ausgabe von graph. Daten steht ein graph. 
SichtgerSt VT-11 (Bildwiederholtyp) mit Lichtgriffel zur Verfiigung. Als 
Betriebssystem wird das Echtzeitbetriebssystem RT-1 1 und als Kom- 
munikationssoftware DECnet/RT von Digital Equipment benutzt; 

Das System wird demnSchst nach Installation einer schnellen Datenver- 
bindung Zugriff auf einen Arrayprocessor ^haben, der bereits an einem 
anderen Projekt eingesetzt wird. Ein Arrayprocessor ist ein schneller 
Spezialrechner, der insbesondere Vektor- und Matrixoperationen sehr' 
schnell durchfilhren kann und ftir den Durchsatz grofier Datenmengen ' 
geeignet ist. Durch Ver legen eines gro(3en Teiles der Auswertepro- 
gramme, insbesondere der rechenintensiven Programmsegmente auf den 
Arrayprocessor konnte die Auswertezeit drastisch (um den Faktor 10-50) 
verkiirzt werden. 


i* ' 

1 . 2 Ein- und Ausgabe graphischer Daten 

Die Eingabe graphischer Daten wird mit einem Kurvendigitalisierer ■ 
(Graf/Pen GP 3 der Fa. SAC), der eine Auflosung von 4000 x 4000 Punkten 
hat, durchgefilhrt. Der Kurvendigitalisierer ist an den Rechner iiber eine 
Koppelelektronik und ein 16-Bit Para lie linterface angesch lessen. 

■^AP 120 B von Floating Point Systems 
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Die graph. Dateneingabe wird benutzt fiir die Erfassung der Interferenz- 
streifenveriaufe bei der halbautomatischen Interferogrammauswertung, 
wobei die Streifen der auf das Tablett des Digitalisierers projizierten 
Interferogramme von einem Auswerter von Hand nachgefahren werden, 
sowie zur Eingabe von Randkonturen der auszuwertenden Objekte. 

Die Eingabe der Daten erfolgt mit einem Datenerfassungsprogramrh 
mit dem im interaktiven Betrieb Funktionen wie Eingeben von Streifen- 
koordinaten und Streifenbezeichnungen, Darstellen der eingegebenen 

Daten auf dem graph. SichtgerSLt in Form von Polygonztigen, Koordinaten- 

i 

transformationen, Durchfiihren einer R.edundanzreduktion der Streifen- 
koordinaten, Loschen von Datenpunkten und Polygonztigen, sowie Ab- 
speichern der Datenstrukturen und Darstellen bereits erfafiter Daten-. 
strukturen auf dem Bildschirm moglich sind. (s. Programm FAW2 , 
Programmbibliothek FAWLIB sowie Programmbeschreibung FAWTXT 
im Tea II) ■ 


1. 3 Ein- und Ausgabe von Grauwertbildern 

Die Eingabe der Interferogramme und Hologramme in den Rechner ge- 
schieht mit Hilfe eines Video-Bilddigitalisierers, der das Videosignal, 
einer hande Istib lichen Videokamera verarbeitet. Es stehen zwei Video- 
kameras mit verschiedenen Aufnahmerohren zur Verfiigung. Eine ist ^ 
mit einer 1" Plumbicon-Rohre und die andere mit einer empfindlicheren 
2/3" Newvicon-Rohre ausgertistet; diese wird zur Aufnahme der relativ 
lichtschwachen Hologrammrekonstruktionen verwendet. Der Bilddigi-'- 
talisierer besteht aus einem programmierbaren Transientrekorder 
(Biomation Model S 8100), der iiber ein Steuerinterface an den Rechner 
angesch lessen ist. Dieses Steuerinterface besteht aus einem Bussystem, 





an das die einzelnen Steuerelemente wie Synchronisationsteil, TV-Zeilen- 
zS-hler und Relaissteuereinheit fiir die Filmprojektorsteuerung angesch les- 
sen Sind. Das Steuerinterface ist modular aufgebaut und konnte leicht 
durch zusatzliche Einheiten (wie z,B. Kamerapositionierung usw, ) erwei- 
tert werden. Eine Skizze des TV -Steuerinterface zeigt Abb, 2. 


Das Video-Signal der Fernsehkamera wird mit dem Transientrekorder 
mit vorgegebener Frequenz (bis zu 100 MHz) mit einer Auflosung von, 

8 Bit digitalisiert und dann vom Rechner per DMA -Transfer ubernommen. 
Der Transientrekorder wird dabei mit dem Synchronisationsimpuls der- 
Fernsehzeile gestartet, bei der die Digitalisierung beginnen soil. Je nach 
gewShlter Abtastfrequenz werden pro Transientrekorderburst zwischen 
1 und 10 Fernsehzeilen erfaCt. Die Auflosung betrSgt in der Zeilenrich- 
tung bis zu 1600 Punkte bei max. ca, 600 Zeilen, wobei allerdings mit 
den zur Zeit verwendeten Videokameras eine Zeilenauflosung von iiber 
1000 Punkten pro Zeile nicht erreicht werden kann. 

Da nach dem Abtasten einer Zeilenfolge bis zum Abtasten der nSchsten- 
Zeilenfolge wegen des Zeilensprungverfahrens 20 ms bzw. 40 -ms gewartet 
werden milssen, bleibt genug Zeit fhr den Datentransfer in den-Rechner 
und fUr einfache Rechenoperationen (wie z.B. Grauwerttransformationen, 
Signa Imitte lung) . 


Fiir die von uns norma lerweise benutzten Auflosungen von 512 Spalten und 
256 Zeilen bzw. 512 Spalten und 512 Zeilen ergeben sich Abtastzeiten 
von ca. 5 sec. bzw. ca. 10 sec. Diese Zeiten kdnnten auf etwa ein Drittel 
reduziert werden, wenn in einem TV -Frame drei Zeilenfolgen abgetastet 
wiirden, was wegen der Datentransferzeit in den Rechner von ca. 5 ms 
moglich ist. 
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Abb. 2: Uber sichtsskizze TV - Interface 










Die Unterprogramme zur Steuerung des Bilddigitalisierers und des Pro- 
jektors, sowie Unterprogramme zur schnellen Bi Iddatenmanipulation 
wie Addition, Grauwerttransformation und zum Datentransfer in den 
Virtual-Memory-Bereich (Speicherbereich von 56 - 248 KB der vom 
Betriebssystem norma lerweise nicht benutzt wird) sind in der Programm- 
bibliothek TVLIB zusammengestellt (s. Teil II), Diese Programme sind 
iiberwiegend im Assemblercode geschrieben, um Verarbeitungsgeschwin- 
digkeiten moglichst klein zu halten und einige hardwarenahe Programme 
einfach zu realisieren. AUe Assemblerunterprogramme sind jedoch von 
FORTRAN-Programmen aufrufbar, so daft die Realisierung komplexerer 
Bildauswerteprogramme in der einfacheren Programmiersprache 
FORTRAN erfolgen kann, ohne da(3 Zeitnachteile in Kauf genommen 
werden miissen. 

Die Ausgabe von Grauwert- oder Binarbildern ist im wesentlichen nur 
in der Entwicklungsphase des Systems notig, um Auswirkungen bestimmter 
Bildverarbeitungsalgorithmen zu untersuchen und zu Dokumentationszwecken. 

Hierzu wird auch das graph, Sichtgerdt VT-11 verwendet, Auf dem Sicht- 
gerS-t kdnnen 1024 x 1024 Punkte adressiert werden, es konnen jedoch 
maximal ca, 2000 Punkte gleichzeitig dargestellt werden, wobei durch 
die Nachleuchtdauer des Schirms sich bei der Anzeige von Binarbildern 
der Eindruck eines stehenden Bildes ergibt. Die Ausgabe von Grauwert- 
bildern geschieht durch je nach Helligkeit mehrfaches Anzeigen der Bild- 
punkte auf dem Schirm, und durch Abfotografieren des Bildschirms mit 
z.B, einer Polaroidkamera, Damit die Anisgabezeit nicht zu groG wird 
(bei der Intensitdt 256 rniiGte der Punkt ja- 256 mal angezeigt werden), 
werden die 8-Bit,Grauwerte linear auf 4-Bit-Werte abgebildet gemaGder 
Transformation 
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wobei g(x,y) der Grauwert an der Stelle (x,y) ist. Die Werte a und g^, 
die man z.B. aus dem Grauwerthistogramm (der prozentualen Verteilung 
der einzelnen Grauwerte) gewinnen kann, konnen dabei zur Ausnutzung 
des Dynamikbereiches vorgegeben werden. 

Da die Punkte des Sichtgerates in einem quadratischen Raster angeordnet 
Sind, die Bildelemente der digitalisierten Bilder aber je nach Wahl von 
Abtastfrequenz und Zeilenabstand nichtquadratisch sind, kann bei der 
Grauwertausgabe iiber das VT-11 Slchtgerat eine Bildstauchung auftreten, 
die jedoch bei der Analyse von Interferogrammen nicht stort. 

Zur Ausgabe von Grauwert- und Binarbildern, sowie zur Berechnung 
von Grauwerthistogrammen dienen die Programme TVDS6, TVDS7A 
und TVHIS, sowie Unterprogramme in der Programmbibliothek TVLIB 
(s. Teil II;. 
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Abb. 3: Beispiele auszuwertender Interferogramme 

a, b) digitalisierte Interferogramme (AuflSsung 512 x 256 Punkte) 

c, d) Verlauf der Intensitat in Zeilenrichtung am Beispiel einiger Zeilen. 
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2 . Interferogrammauswertung 

Ein Interferogramm 1st eine Darstellung des untersuchten Objektes 
durch eine Schar von Linien konstanter Lichtweglange, die durch Inter- 
ferenz zweier Lichtwellen entstehen. Bei der Doppelbelichtungsholo- 
grafie werden zwei Lichtwellen, die zwei Zust^den des Testobjektes 
entsprechen, auf einer Ho logrammp latte gespeichert, das Interferenz- 
muster ergibt sich bei der Rekonstruktion. Es gibt im wesent lichen 
zwei Verfahren zur Auswertung von holografischen Interferogrammen, 
namlich das statische und das dynamische Verfahren [9]. Beim stati- 
schen Verfahren wird das virtuelle rekonstruierte Bild des Objektes 
durch die Hologrammp latte von einem festen Beobachtungspunkt aus 
betrachtet und die relative Interferenzordnung an jedem Punkt durch ' 
AuszShlen der Streifen bestimmt, Beim dynamischen Verfahren wird 
die Anderung der Interferenzstreifenordnung an einem festen Punkt 
im virtuellen rekonstruierten Bild des Objektes bei Ver^nderung des 
Beobachtungspunktes ermittelt. Das dynamische Verfahren ist jedoch 
nur dann anwendbar, wenn die Interferenzstreifen nicht auf dem Objekt 
lokalisiert sind. Von uns wird zur Auswertung der holografischen 
Interferogramme das statische Verfahren angewendet, 

Im Falle der StrSmungsinterferogramme handelt es sich bei den Inter- 
ferenzstreifen um Linien gleicher mittlerer Dichte des Stromungsme- 
diums, da der Brechungsindex von Gasen direkt mit der Dichte zusam- 
menha.ngt, wShrend die Interferenzstreifen in den Reifentesthologram- 
men Linien gleicher Hohe, d. h. Linien gleicher Objektauslenkung sind. 
Bei diesen Interferogrammen handelt es sich um Zweistrahlinterfero- 
gramme, was zu einer sinusformigen IntensitStsmodulation fiihrt, die 
im Falle der Reifentesthologramme noch von der Objektintensitat iiber- 
lagert sind. 


Ziel der Auswertung dieser Interferogramme ist die Rekonstruktion 
der Objektfunktion, d. h. der Dichtefunktion beiden Stromungsinter- 
ferogrammen und der Verformungsfunktionen bei den Reifentestholo- 
grammen. Bei diesen sind auBerdem noch bestimmte Sirukturen 
in den Interferenzmustern zu erkennen und zu lokalisieren, durch die 
auf das Vorliegen bestimmter Reifenfertigungsfehler gesch lessen wer- 
den kann. Fiir beide Arten der Auswertung ist das Erfassen der Inter- 
ferenzstreifenver laufe erforderlich. Fiir die Rekonstruktion der Dichte- 
funktion Oder Verformungsfunktion ist auBerdem die relative Numerie- 
rung der Streifen, d. h. die Ermittlung der relativen Lichtweglange 
eines jeden Streifens erforderlich, um daraus auf die absolute Gr5Be 
der Dichte bzw. der Verformung umrechnen zu konnen. 

Die Erfassung der Interferenzstreifen geschieht bei der sog. halbauto- 
matischen Auswertung mit Hilfe eines Kurvendigitalisierers durch 
Nachfahren der Streifenver laufe von Hand (s. a. Kap. 1.2). Bei der 
vollautomatischen Auswertung muB der Streifenver lauf der Interferenz- 
streifen aus dem digitalisierten Grauwertbild, eventuell nach Anwendung 
geeigneter Bildvorverarbeitungsschritte, vomRechner extrahiert war- 
den. 

Nach diesem Schritt liegen die Interferogramme, d.h. ihre Interferenz- 
ordnungsfunktionen in Form einer graphischen Datenstruktur imRechner 
vor, bei der jeder Streifen durch zwei Polygonziige, die den Ver lauf 
seiner beiden Flanken darstellen, reprasentiert ist. Diese Darstellung 
wird nun auf das Vorhandensein von Streifenunterbrechungen und Strei- 
fenfehlverbindungen, die durch Kontrastschwachen oder zu hohe lokale 
Streifendichten in den Hologrammen entstehen konnen, gepriift und evtl. 
verbessert. 
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Daran schliefit sich die Streifennumerierungsphase an, die so ausgelegt 
ist, dafi hier evtl. in den vorhergehenden Schritten nicht erkannte Strei- 
fenfehlverbindungen erkannt und aufgeldst werden. 

Interferenzlinien haben die gleichen Eigenschaften wie Hohenlinien auf 
einer Landkarte. Aus dieser Eigenschaft resultieren einige wichtige 
Regeln fiir die Streifennumerierung, weitere Regeln ergeben sich aus 
der Natur des Objektes. 

Aus dem Vergleich der Interferenzstreifen mit HShenlinien geht schon 
hervor, dafi die Datenstruktur , in der die Interferogramme nach diesem 
Schritt vorliegen, fiir die weitere numerisch-mathematische Auswertung 
nicht geeignet ist. Bei den Stromungsinterferogrammen sollen zum Bei- 
spiel Zeitreihen der Dichte aus einer Reihe von Interferogrammen an 
vorgegebenen Orten des Mefifeldes gebildet werden, um daraus Kreuz- 
korrelationen und PhasenplSne zu berochnen. Im Falle von turbu lenten 
Stromungen ist aufierdem eine Mittelwertbildung iiber mehrere Inter- 
ferogramme notig, um statistische Schwankungen zu eliminieren und 
Vergleiche mit theoretischen Ergebnissen zu erzielen. Im Falle der 
Reifentesthologramme ist u. a . eine Untersuchung von Gradienten und 
Kriimmungen der Verformungsfunktion, sowie die Ermittlung einer 
Normverformungsfiache und die Bildung von Differenzinterferogrammen 
vorgesehen. 

A lie diese Operationen erfordern, dafi die Interferenzordnungsfunktionen 
in einer Form F = f(x,y) vorliegen. Zur Durchfiihrung dieser Transfor- 
mation wurde ein Verfahren der gewichteten zweidimensionalen Polynom- 
approximation mit nachfolgender Splineinterpolation benutzt, wodurch die 
Interferenzordnungsfunktionen als Satz von Polynomkoeffizienten an geeig- 
net gewahlten Stiitzstellen dargestellt werden. 
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In den folgenden Abschnitten werden die einzelnen Verfahrensschritte 
wie Bildvorverarbeitung, Interferenzstreifenextraktion, Streifennume- 
rierung und Polynomapproximation genauer dargestellt. Es schliefSt 
sich ferner ein Kapitel an, welches die Erkennung spezieller Strukturen 
in den Reifentesthologrammen beschreibt, sowie ein Kapitel iiber die 
Entzerrung von Objektgeometrien. 


2. 1 Bildvorverarbeitung 

Nach der Digitalisierung der Vorlagen ist es in vielen Fallen erforder- 
lich vor der Extraktion der Interferenzstreifen eine Bildvorverarbeitung 
vorzunehmen. Dazu zShlen raumliche und zeitliche Mittelungsverfahren, 
Gradientenbildung, Subtraktion der mittleren Intensitatsverteilung, An- 
wendung von zweidimensionalen Filterprozeduren. 


2.1.1. Raumliche und zeitliche Mittelung 

Die digitalisierten Interferogramme enthalten einen gewissen Rauschan- 
teil, der einerseits bei der Hologrammerstellung und -reproduktion 
durch Speckle-Bildung und Streuung, andererseits durch die elektronischen 
Komponenten bei der Interferogrammaufnahme und -digitalisierung entsteht. 

Der Rauschanteil , der durch das Kamera-Rauschen und die Digitalisierung 
entsteht, ist norma Ivertei It und hat ein Signal-Rausch-Verhaitnis von ca. 

36 dB (bei 0,4 V Eingangsempfindlichkeit). Daher kann eine Verbesserung 
des Signals insbesondere bei sehr kcntrastarmen Hologrammen durch 
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Mittelwertbildung iiber mehrere zeitlich nacheinander aufgenommene 
Zeilensignale (sog. Bildstapelung) erreicht warden. Bei den hologra- 
fischen Interferogrammen konnte das Signal-Rausch-VerhSltnis durch 
Signalmittelung iiber 20 Bilder um ca. 13 dB verbessert warden. Das 
orts- und zeitfeste Rausch en (Speckles) wird durch Anwendung einer 
raumlichen Mittelung, wobei der Grauwert an jedem Punkt des Bildes 
durch das arithmetische Mittel seiner 8 Nachbarpunkte ersetzt wird, 
vermindert. Ein Beispiel zur zeit lichen und rSumlichen Mittelung 
zeigt Abb. 5 (s. a. Programme TVMN6, TVMN7A und SM08 in Teil II). 
Bei Anwendung von raumlichen Mittelungsverfahren wird allerdings die 
Ortsauflosung vermindert, was insbesondere bei hohen Streifendichten 
beachtet warden muB. 


2.1.2. Subtraktion der UntergrundintensitSt 

Wenn Interferenzstreifen in den Hologrammen mit geringem Kontrast 
an Orten verlaufen, wo die Untergrundintensitlt einen grofieri raumlichen 
Gradienten hat, kann es vorkommen, daB die Streifen mit dem im Ab- 
schnitt 2. 2 diskutierten Extraktionsverfahren nicht mehr erkannt warden, 
obwohl das Verfahren im Prinzip eine iiber lagerte Untergrundintensitat 
beherrscht. Um die Resultate mit dem Extraktionsprogramm zu ver- 
bessern, wurde ein Vorverarbeitungsprogramm entwickelt, das eine 
FlMche mittlerer Intensitat durch Mittelwertbildung in einem Raster von 
32 X 32 Punkten und linearer Interpolation bildet und vom Original sub- 
trahiert. Ein Beispiel zur Subtraktion der mittleren Grauwertfiache 
zeigt Abb, 6 (s.a. Programme SUBMF, SUBMF3 im Teil II). 


n*3 


Abb. 5; Beispiel zur raumlichen und zeitUchen Mitte’ung 

a) ungemittelte Zeile aus einem Hologramm (512 Bildpunkte) 

b) die gleiche Zeile nach Mittelung iiber 5 zeitlich hinterein- 
ander aufgenommene Signale 

c) wie b) aber iiber 20 Signale gemittelt 

d) die raumlich (iiber 8 Xachbarpunkte) gemittelte Zeile a 
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Abb. 6: 


j 


Hologramm Nr. 04 Zeile 350 


d 


e 


f 


Hologramm Nr. 18 Zeile 100 






Zeilenprofile zur Demonstration der Subtraktion der Unter- 
grundintensitat (Programm SUBMF 3) 

a, d) Zeile aus einem Hologramm 

b, e) Zeile aus der berechneten Fl&che mittlerer Intensitat 

c, f) Hologrammzeile nach Subtraktion der mittleren Inten- 

sitatsfiache 



Falls man zur Aufnahme und Reproduktion der Hole gramme ein Ver- 
fahren benutzen wiirde, welches zwei oder mehr unterschiedliche Refe- 
renzstrahlen verwendet, kbnnte man die Unterdriickung der Untergrund- 
intensita.1 und der Streuerscheinungen an Staubteilchen durch Differenz- 
bildung zweier Signale erreichen, die man durch eine Phasenverschie- 
bung der Referenzstrahlen bei der Rekonstruktion des Hologramms er- 
halt. 


Bei den Hologrammen der Fa, Continental sind solche Verfahren aller- 
dings nicht ohne weiteres anwendbar, da dazu grundlegende Anderungen 
an der Holografieanlage durchgefiihrt werden mufJten, 


2.1,3. Storbeseitigung durch zweidimensionale Filterung 

Eine weitere in den Reifentesthologrammen auftretende StSrung ist das 
durch die Heizbalgabdrticke hervorgerufene iiberlagerte Stormuster. Die 
Storung des Interferenzmusters durch die Heizbalgabdriicke ist teilweise 
recht groB (s. Abb. 7) und fiihrt einerseits dazu, daB schon bei der 
Streifenextraktion Fehler auftreten, andererseits wird die nachfolgende 
Reifenfehleranalyse erschwert, da die durch das Storinterferenzmuster 
erzeugten Strukturen in GroBe und Form nicht von solchen Strukturen 
zu unterscheiden sind, die auf Reifenfehler zuriickzuftihren sind. Daher 
ist es erforder lich, dieses Stormuster vor der Streifenextraktion zu eli- 
minieren. Bei den Heizbalgabdriicken handelt es sich um schrag zur Reifen- 
laufrichtung parallel zueinander verlaufenden Strukturen, die durch die 
Art der Hologrammaufnahme und Rekonstruktion im digitalisierten Bild 
nicht genau aber annahernd als paralleles gerades Streifensystem er- 
scheinen, das dem eigentlichen Streifenmuster iiberlagert ist. 
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Auf Grund der regelmafligen Struktur des Stormusters stellt die zwei- 
dimensionale Fourieranalyse eine Moglichkeit dar, das StSrmuster zu 
eliminieren, 

Wie aus der Theorie der Fourieroptik bekannt ist, entspricht die im 
Fernfeld beobachtete Beugung von parallelem Licht an einer zweidimen- 
sionalen Struktur einer Fouriertransformation dieser Struktur [12]. 

Optisch kann man das Fernfeld durch Abbilden mit einer Linse in der 
Brennebene der Linse beobachten. Beugt man z.B. paralleles Licht an 
einem Gitter mit dem Linienabstand a, so erhalt man als Beugungs- 
bild Punkte, die mit einem Abstand proportional zu 1/a auf einer Geraden 
liegen, die senkrecht zu den Gitterlinien orientiert ist. Die Intensitat 
der Punkte entspricht der Transmissionskurve des Gitters und stellt 
bei einer rechteckformigen Transmission eine sin(x)/x-Verteilung dar. 

Hat man ein Bild in der Form f(x, y), wobei f(x, y) dem Grauwert an 
der Stelle (x, y) entspricht, so ist die Fouriertransformation F(u, v) 
definiert als 

00 

F(u, v) = JJ f(x, y) exp (- i 2 7T (ux + vy)) dxdy 
- 00 

und die Rticktrans formation 

CO 

f(x, y) = JJ F(u, v) exp (i2 7r(ux + vy)) dudv 
- 00 

Bei der zweidimensionalen Fourieranalyse werden die Koeffizienten einer 
Entwicklung nach zweidimensionalen Sinus- bzw. Kosinusfunktionen be- 
rechnet, die den relativen Anteil von Grundfunktionen der Form 
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exp(i 2 71 (ux + vy)) am betreffenden Bild darstellen. Man nennt daher 
F(u, v) das Frequenzspektrum des Bildes f(x, y). Hohe Raumfrequenzen 
im Bild erzeugen Spektralkomponenten, die an Stellen mit grofJem u bzw. 

V in der Spektral- bzw. Fourierebene liegen. Durch VerSndern bestimmter 
Spektralkomponenten in der Fourierebene konnen die Anteile, die bestimm- 
ten Strukturen in der Bildebene entsprechen, herausgefiltert werden. 

Abb. 7a zeigt ein digitalisiertes Reifentesthologramm, das ein starkes 
Stbrmuster entha.lt. In Abb. 7b ist das Amplitudenspektrum dargestellt. 
Hier sind die Anteile des Stormusters deutlich als he lie Flecken zu er- 
kennen, die auf einer Geraden senkrecht zum Stormuster liegen. Die 
relativ starken Frequenzkomponenten auf der u-Achse sind zum groBen 
Teil auf das Erweitern der 480 Zeilen des digitalisierten Bildes auf 512 
Zeilen zuriickzuftihren, was auf Grund der Anwendung von Fast-Fourier- 
Transform-Algorithmen [13] erforderlich ist und zunachst willkurlich 
durch Nullsetzen der zusatzlichen Zeilen geschah. Dadurch ergibt sich 
naturlich in y-Richtung eine Unstetigkeit an der Zeile 480, was zu starken 
Spektralkomponenten in dieser Richtung fiihrt. Durch Multiplikation der 
Spektralebene mit einer geeignet gewahlten Filterfunktion und Riicktrans- 
formation des modifizierten Spektrums erhait man das gefilterte Bild. 

Die Filterfunktion wird zunachst durch Nullsetzen der entsprechenden 
Spektralkomponenten realisiert. Hierbei werden auch gleich die Spektral- 
komponenten, die hohen Raumfrequenzen im Bild entsprechend auf Null 
gesetzt, um eine BildglSttung zu erreichen, die sonst durch Anwendung 
von Mittelungsoperatoren in der Bildebene erfolgen miifite. 

Abb. 7d zeigt das Hologramm nach Riicktransformation aus dem modifi- 
zierten Spektrum der Abb. 7c. Man erhS.lt das urspriingliche Hologramm 
zuriick, es ist jedoch das Stormuster weitestgehend beseitigt. 

Zur Durchfiihrung der Fourierfilterung dienen die Programme ZWDFFT, 
ZWDCLR, ZWDRFT, ZWDAMP und TRANSP (s. Teil II). 
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Abb. 7: Storbeseitigung durch Filterung in der Fourierebene 

(Hologramm Nr. A4 8007228) 

a) Hologramm mil iiberlagertem Storinterferenzmuster 
durch Heizbalgabdriicke 

b) Spektrum des Hologramms 

c) Gefiltertes Spektrum 

d) Rticktransformation des gefilterten Spektrums 





c, o 



Es ist neben dem Abschneiden der hohen Raumfrequenzen auch moglich, 
die niedrigen Raumfrequenzen, die ja im we sent lichen die Untergrund- 
helligkeit in den Hologrammen darstellen, zu beeinflussen. 

Erste Ergebnisse mit dem Ausfiltern der Grundfrequenzen haben jedoch 
gezeigt, dafi es offenbar problematisch ist, diese Spektralkomponenten 
einfach auf Null zu setzen, da durch die dann vorhandene grofie Unstetig- 
keit der in diesem Bereich sehr gro(3en Werte der Spektralkomponenten 
nicht ganz befriedigende Resultate bei der Riicktransformation auftraten. 

Solche zweidimensionalen Filterverfahren erfordern recht hohe Rechen- 
zeiten, wenn sie mit dem ProzefSrechner durchgefiihrt werden (ca. 20 min. 
fiir eine Hin- und Riicktransformation), Nach Inbetriebnahme des Array- 
Processors kann diese Zeit jedoch auf einige Sekunden reduziert werden, 
wodurch ein Einsatz solcher Filterverfahren in der Serienauswertung 
durchaus inFrage kommt. 


2, 2 Extraktion der Interferenzstreifen 


A 


Die Extraktion der Streifen aus der digitalisierten Grauwertdarstellung 
geschieht in zwei Schritten: zuerst miissen aus der lokalen Grauwert- 
verteilung die Objekte, die als Interferenzstreifen interpretiert werden 
konnen, durch eine entsprechende Kontrastanhebung hervorgehoben wer- 
den, anschliefiend miissen die Orte der Streifen, die sowohl aus den 
IntensitMtsextrema wie aus den Flanken abgeleitet werden kbnnen, zu 
zusammenhangenden Linien verbunden werden. Bei ortlich sehr ausge- 
dehnten Streifen sind die Extrema u.U. nicht genau zu bestimmen, so dafi 
die Koordinatenfolgen der ansteigenden und abfallenden Flanken die 
Streifen genauer definieren [10]. 


3 ^ 
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2.2,1 Grauwerttransformation 

Die iiblicherweise in der digitalen Bildverarbeitung benutzten lokalen 
Operatoren zur Kantenverscharfung und Linienerkennung wie Gradient-, 
Laplace- oder Sobeloperatoren, sowie andere komplizierterer Art [14-16], 
Sind im Falle der Interferenzauswertung nicht geeignet, da sie nur auf 
scharfe Grauwertstufen ansprechen und bei den in der Regel kosinus- 
quadratformig modulierten Interferogrammen nicht befriedigend arbeiten. 

Fischer et al. [17] benutzten zur zeilenweisen Extremabestimmung der 
Interferenzstreifen (nach einer Rauschbefreiung) abschnittweise ange- 
paf3te Sinusfunktionen, deren Extrema dann als Streifenorte genommen 
werden. Ein anderes Verfahren, das jedoch parallels Streifen mit kon- 
stantem Abstand, die aufierdem senkrecht zur abgetasteten Zeile ver- 
laufen, voraussetzt, ist in [18] angegeben, Hier wird eine Fouriertrans- 
formation der Intensitat in Zeilenrichtung durchgefiihrt, um die Ortsfre- 
quenz und Phasenlage des Streifenfeldes zu bestimmen, Als Mindest - 
streifenzahl in der Zeile werden dort zehn Streifen angegeben, damit 
ein eindeutiger Peak im Frequenzspektrum entsteht. 

Diese Verfahren sind fiir unsere Auswertung nicht amvendbar, da sie 
im wesentlichen gleichformig gestaffelte Interferenzstreifen, die vor- 
nehmlich senkrecht zur Abtastrichtung verlaufen, voraussetzen. Aus 
diesen Griinden wurde ein von uns zunM,chst erprobtes Verfahren der zei- 
lenweisen Extremabestimmung, das mit einer variablen Schwelle zur 
Rauschunterdriickung arbeitete (s. Zwischenbericht III [1]) modifiziert 
und weiterentwickelt. 

Bei diesem Verfahren werden zundchst solche Extrema in der Zeile ge- 
sucht, die zum zuletzt gefundenen Extremum eine bestimmte Grauwert- 
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differenz aufweisen, die groBer ist als eine Schranke, die sich aus der 
mittleren R aus champ litude ableitet. 

Zur Anpassung an die ortsabh&ngige Kontrastverteilung in den Hologrammen 
wurde dort statt einer festen Schranke eine ortsabhM,ngige Schranke be- 
nutzt (s.a. Programme STREX3 und STREX4 im Teil II). Die Extrema 
werden nun benutzt, um eine treppenformige Schwellwertfunktion zu 
definieren, die jeweils zwischen benachbarten Extrema die mittlere Am- 
plitude annimmt. Am Zeilenanfang bzw. am Zeilenende wird die Trigger- 
schwelle der ersten beiden bzw, der letzten beiden Extrema benutzt. Der 
Schwellwertvergleich der Intensitatswerte mit dieser Treppenfunktion 
liefert dann eine Binardarstellung der betrachteten Zeile (siehe Abb, 8). 

Abb. 9 zeigt einen Ablaufplan des Verfahrens. Diese Art der Kontrast- 
anhebung entspricht einer Flankentriggerung der Streifen nach einer ent- 
sprechenden Elimination der UntergrundintensitS.t, wie sie auf elektroni- 



Abb. 8: 


Zur Grauwerttransformation 

(oben: Grauwertzeile 

unten: daraus abgeleitete BinM.rzeile) 



nt f? c 1 1 s I. o n 



Ablaufplan zur zeilenweisen Extremasuche 
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schem Wege [10], oder durch Anwendung des Programms SUBMF3 
(Kap. 2,1,2) mit einer fasten Schwelle durchgefiihrt werden konnte. 

Die Ergebnisse der zuletzt b esprochenen Kontrastanhebung sind jedoch 
im allgemeinen wesentlich besser, als wenn lediglich die Untergrund- 
intensitat mit dem Programm SUBMF3 eliminiert wind. Die besten 
Resultate wurden erzielt, wenn nach Anwendung des Programms SUBMF3 
die Kontrastanhebung angewendet wurde. Die Abb, 10 zeigt eine An- 
wendung des Programms STREX3 auf zwei Hologramme, Man sieht, dafi 
das Verfahren selbst bei sehr stark variierender Untergrundhelligkeit 
und schwachem Kontrast eine gute Trennung der Interferenzstreifen 
vom Untergrund ergibt, Trotz sehr guter Rauschunterdrtickung sind 
auch kleine Strukturen noch voll im Interferogramm enthalten. 


2, 2, 2 Polygonzugdarstellung 

Das nach Kap, 2,2,1 in eine Binarform gebrachte Interferogramm wird nun 
in eine Polygonzugdarstellung umgeformt, die die Koordinatenfolgen der 
StreifenrSnder darsteUt, Dadurch wird eine Datenreduktion von ca, 30 : 1 
bis 100 ; 1 erreicht, Es werden jeweils zwei aufeinanderfolgende Zeilen 
des BinMrbildes miteinander vergUchen und durch die Priifung von Fort- 
setzungsbedingungen [19] entschieden, wo ein Streifen aus der alten Zeile 
in der neuen fortgesetzt wird, ob ein Streifen verzweigt, neue Streifen ent- 
stehen oder alte verschwinden, Durch das sequentielle zeilenweise Abar- 
beiten des Bildes ist ohne weitere Prtlfung sichergestellt, dafi alle Objekte 
im Bin&rbild erkannt wurden, was bei einer Randverfolgung der einzelnen 
Objekte nicht einfach ist, AuBerdem braucht auf jede Bildzeile nur einmal 
zugegriffen zu werden, wodurch die Verarbeitung direkt im AnschluB an 
die Digitalisierung durchgefiihrt werden kann und eine Speicherung des ge- 
samten Bildes nicht erforderlich ist. 




Abb. 10: Beispiele zur Grauwerttransformation 

a) Grauwertdarstellung des digitalisierten Hologramms 
Nr. 16 auf dem graphischen Sichtgerat 

b) Binardarstellung des Hologramms nach Bildverbesserung 
und Anwendung des Programms STREX 3 

c) wie a) Hologramm Nr. 05 

d) wie b) Hologramm Nr. 05 
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Die Streifenkoordinaten der linken und rechten Objektr finder werden ent- 
sprechend der Streifenzugehorigkeit unter einer Adresse, die dem je- 
weiligen Streifenrand zugeteilt wird, abgespeichert. Damit nicht zu viele 
redundante Koordinaten gespeichert werden, wird eine Redunanzreduktion 
nach [20] durchgefiihrt, Mit dieser Reduktion wird die quasi kontinuier- 
liche Koordinatenfolge durch Polygonziige mit minimaler Anzahl von 
Koordinaten ersetzt. Die Genauigkeit, mit der die ursprunglichen Koor- 
dinatenfolgen durch die Polygonziige dargestellt werden, wird durch eine 
Toleranzschranke bestimmt. 

Durch die Verfolgung von linken und rechten Streifenrfindern wird die 
Streifendichte in den Interferogrammen verdoppelt, was sich aber ins- 
besondere bei breiten Streifensystemen giinstig auf die weitere Verar- 
beitung (Berechnung der Objektfunktion) auswirkt. Wiirde man zu dem 
Verlauf der Flanken auch noch den Verlauf der Extrema mit erfassen, 
so konnte das Interferogramm durch eine vierfache Liniendichte noch 
genauer dargestellt werden. Wenn wir hier von linken und rechten Streifen- 
rfindern sprechen, so beziehen wir uns dabei auf die abgetastete Zeile. 

Bei geschlossenen Streifen oder solchen, die mehrmals von einer Zeile 
geschnitten werden, geht der linke Streifenrand in den rechten fiber und 
umgekehrt, 

Ein Problem, das bei der Streifenverfolgung auftritt, ist die Erkennung 
und Behandlung des Bildrandes, Endet ein Objekt des Binfirbildes am 
Rand, so dfirfen die beiden davon betroffenen R finder nicht miteinander 
verbunden werden, da sie Linien unterschiedlicher Ordnung sein konnen. 
Das gleiche Verbindungsproblem ergibt sich, wenn Objekte im Bild enden. 
Stellt das Objekt einen geschlossenen Interferenzstreifen dar, so sind die 
beiden Rander zu verbinden. Endet ein Objekt, welches einen sehr schma- 
len Streifen darstellt, weil z.B. die Auflosungsgrenze unterschritten wird. 


so diirfen die Render nicht verbunden werden. Dieses Problem stellt 
eine prinzipielle Schwierigkeit dar, die nicht allein aus der Analyse 
der lokalen Grauwertverteilung in der Umgebung des Ortes geklSrt 
werden kann, sondern es miissen die globalen Informationen iiber den 
Verlauf des Streifenfeldes in die Entscheidung iiber den Zusammen- 
hang solcher Gebiete mit eingehen. Fiir den Fall der Bildberandung 
wird das Problem durch Vorgabe der Geometrie des Objektrandes und 
einem entsprechendem Randtest bei der Streifenverfolgung gelost. Zur 
Umformung des Interferogramms in die Polygonzugdarstellung wird 
das Programm STR2 und zur Vorgabe der Geometrie des Objektrandes 
das Programm IVIASKOl zusammen mit dem Programm FAW2 benutzt 
(s.a. Programmlisten der Programme im Teil II). 

Beispiele fiir die Extraktion der Interferenzstreifen aus dem Grauwert- 
bild nach der KontrastverstSrkung zeigen die Abb. 12 - 17. 


2. 2. 3 Behandlung von Streifenfehlverbindungen 

Bei der Extraktion der Interferenzstreifen aus dem Hologramm und der 
Darstellung als Polygonziige erhilt man gelegentlich fehlerhaft verbundene 
Linien, wenn Streifen durch KontrastmSngel im Hologramm, oder durch 
Unterschreiten der Auflosungsgrenze nicht sicher erkannt werden. 

Es ist nun wiinschenswert, den GroiBteil dieser Fehlverbindungen schon vor 
Anwendung der Numerierungsverfahren zu berichtigen. Deshalb haben wir 
das Programm BVSF (s. Teil II) entwickelt, das gezielt nach Stellen sucht, 
wo Interferenzstreifen innerhalb des Bildes enden. Es wird versucht, durch 
VerlSngerung dieser Streifen auf einen ebenfalls vorzeitig endenden Inter- 


ferenzstreifen zu treffen und die Streifenunterbrechung durch Verbindung 
entsprechender F lankenlinien zu beseitigen, Wie an dem Beispiel in Abb. 11 
zu erkennen ist, konnen auf diese Weise die meisten Fehler dieser Art be- 
seitigt werden. Von dem Verfahren nicht erkannte Streifenunterbrechungen 
und Linienfehlverbindungen werden bei der Numerierung des Linrenbildes 
erkannt und beseitigt. 

Fine weitere Methode zur Lokalisation von fehlerhaft extrahierten F lanken- 
linien, das im Wesentlichen die geometrischen Parameter des Streifens, 
bzw. der vom Streifen eingeschlossenen FlM.che benutzt, wird in Kap. 2, 6 
beschrieben. 


2. 3 Streifennumerierung 

Die Numerierung von Interferenzstreifenbildern unterliegt gewissen Regeln, 
die teils allgemein gelten und sich teils aus der Natur des untersuchten 
Objektes ergeben. Diese Regeln sind schon im Zwischenbericht IV [1] 
beschrieben worden, auBerdem wurde dort eine Strategie erklart, mit der 
einfache Interferogramme automatisch zu numerieren sind. 

Hiernach werden zunachst die inneren Linien von Liniensystemen, sowie 
die inneren Linien von am Bildrand endenden Liniensystemen aufgesucht. 
Ausgehend vom Liniensystem mit den meisten Linien werden nun die 
Schnittpunkte der Verbindungsgeraden zwischen den Flachenschwerpunkten 
der inneren Linien der Systeme gebildet, und die Streifen entsprechend 
der Reihenfolge der Schnittpunkte numeriert. Die Linien, die nach Abar- 
beiten der einzelnen Liniensysteme noch nicht numeriert werden konnten, 
werden durch Betrachten der Liniennachbarn numeriert. Diese Strategie 
wurde von uns erfolgreich auf die Numerierung von fehlerfrei extrahierten 
Streifenfeldern angewendet. 
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Um jedoch auch Streifensysteme mit fehlerhaft erkannten Linien nume- 
rieren und berichtigen zu konnen, was mit dem oben beschriebenen Ver- 
fahren zu unvertretbar hohem Programmieraufwand gefiihrt hatte, ent- 
schieden wir uns ftir die Entwicklung einer zweiten Numerierungsstrate- 
gie (s, Programm NUMl und Programmbibliothek NUMLIB in Teil II). 
Eine wesentliche Erkenntnis, die zu dieser Strategic gefiihrt hat, ist die 
Erfahrung, da(3 KontrastschwSchen und hohe Liniendichten, die zur Sto- 
rung des Linienver laufs und des Linienzusammenhangs in den extrahier- 
ten Streifenfeldern fiihren konnen, hSufig nur an wenigen Stellen des 
Bildes (vornehmlich im Randbereich bei den Hologrammen) auftreten, 
wahrend der Linienver lauf auBerhalb dieser Stellen einwandfrei erkannt 
wird. 

Zur Lokalisation dieser Storstellen wird das Linienbild des Interfero- 
gramms durch ein geeignet gewShltes Maschennetz in einzelne Segments 
zerteilt. Ftir jede Masche wird nun eine Zahl (Konfidenz) ermittelt, die 
Aufschlufi gibt iiber die Numerierbarkeit der einzelnen Masche. In die 
Konfidenzzahl geht neben der Anzahl von verschiedenen Linien in der 
Masche auch die Anzahl von Stellen ein, an denen sich die Richtung der 
Liniennumerierung umkehrt. Als Ausgangspunkt fiir die Numerierung 
wird nun die Masche ausgew§.hlt, die die wenigsten Umkehrzentren und 
die meisten Linien besitzt. Hat diese Masche genau zwei Umkehrzentren, 
so ist sichergestellt, dafi es sich um eine fehlerfreie Masche handelt. 

Die Numerierung wird nun bei einer der vier Nachbarmaschen fortge- 
setzt, die von den in Frage kommenden Maschen die hochste Konfidenz 
hat und in der die Numerierung fortgesetzt werden kann, das heifit welche 
mit der ersten Masche mindestens zwei Schnittpunkte verschiedener 
Linien gemeinsam hat. Die Nummern der anderen Nachbarmaschen wer- 
den in einem Stack abgespeichert, aus dem im folgenden die Maschen- 
nummern entnommen werden, welche als nS.chste zu numerieren sind. 
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1st dieser Stack abgearbeitet, aber noch nicht alle Maschen numeriert, 
so werden, von der schon numerierten Masche mit der besten Konfi- 
denz ausgehend, sog, Grofimaschen gebildet, die aus zwei oder mehr 
benachbarten Maschen bestehen und die Numerierung so fortgesetz^t. 

Stofit das Numerierungsverfahren bei einer Masche auf einen Wider - 
spruch, so wird sie zunSchst zurtickgestellt. Nach Abarbeiten alter 
anderen Maschen werden die fehlerhaften Maschen bearbeitet. Es wer- 
den jetzt die Nummern, die am Rand der Masche von der Numerierung 
der Nachbarmaschen her bekannt sind, ilbernommen, sofern sie mit dem 
Linienver lauf in der Masche vertrtglich sind. Alle Streifen, die mit 
zwei verschiedenen Nummern versehen werden miifiten, werden in der 
Masche aufgetrennt und jedes einzelne Segment wird mit der Linien- 
nummer aus der Nachbarmasche versehen. Durch diese Strategie wirken 
sich Storungen im Linienver lauf nicht auf die Numerierbarkeit des 
ganzen Liniensystems aus, sondern allenfalls auf die nSchste Nachbar- 
schaft. 


2. 4 Polynomdarstellung 

Die Durchfiihrung einer Transformation der in Polygonzugform vorlie- 
genden Interferenzordnungsfunktion in eine Darstellung der Form z = f(x,y) 
ist gleichbedeutend mit dem Problem.eine Funktion aus ihren Konturlinien 
(Linien mit f(x,y) = const. ) zuriickzugewinnen. 

Die dabei notwendige Interpolation auf einer durch beliebig angeordnete 
StUtzwerte vorgegebene FlSche ist schwieriger durchzufUhren, als wenn 
die Stiitzwerte der FlSche auf einem regelma(3igen Gitter oder zumindest 
auf parallelen Linien angeordnet wSren. Filr diese regelm^(3igen Stiitz- 
punktverteilungen gibt es auch praktisch akzeptable Verfahren, die ortho- 
gonale Polynome [21-23] oder Splinefunktionen benutzen [24, 25]. 
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Ein von Hayes und Halliday [24] benutztes Verfahren, das beliebige 
Stiltzpunktverteilungen beherrscht, wobei in einem regelma(3igen, tiber 
die Stiitzpunkte gelegten Gitter sog. B -Splines berechnet werden, filhrt 
zu einem unvertretbar hohen Rechenaufwand (s. a. Zwischenbericht II [1]). 

Nach der von uns verwendeten Methode werden die Stiitzpunkte in der 
NShe der Interpolationsstelle mit ihrem Abstand gewichtet und durch 
ein Polynom niedrigen Grades approximiert. Den Wert an der Interpola- 
tionsstelle erhSLlt man nun durch Auswerten dieses Polynoms [26], 

Dieses Verfahren der lokalen Polynomapproximation wurde ftlr die 
Darstellung der Interferogramme gewShlt, weil es relativ kurze Rechen- 
zeiten erfordert und durch Auswerten der Polynomkoeffizientenmatrix 
sofort wichtige Informationen wie Stellen grofier Kriimmung oder Stellen 
groI3er Gradienten zu erhalten sind. AuBerdem ist eine leichte Beseiti- 
gung von ilber lagerten Keilfeldern mbglich. 

Das AieBfeld wird dazu von einem regelmaBigen Gitter mit den Gitter- 
punkten ( A Pj) (i = 1, . , n, j = 1 , . . m) iiberzogen. An den Gitterpunkten 
werden nun zugunsten des Rechenaufwandes Polynome zweiten Grades 
angesetzt: 

2 2-r 

(-) P(x. y) = 2 2 ''rs'''' 

r=0 s=0 


Die Koeffizienten an der Stelle ( ,a., p^) werden mit der Nebenbedin- 

gung bestimmt, daB die Summe der Fehlerquadrate minimal ist (GauB’sche 

Methode der kleinsten Quadrate); 

n 

(=-•) Q = 2 ■ ^r^^ ^ i^^ ^ ^^r 

r=l 
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w(d ) ist eine Gewichtsfunktion, die die Sttitzwerte mit grofiem Ab- 

stand zu ( A fi.) geringer gewichtet. Durch die Verwendung einer 

^ ^ 2 2 
stark abfallenden Funktion fiir w(d ) wie exp( - ct d ) hat man den 

Vorteil, nicht alle Stiitzpunkte zur Berechnung von P.^ 

bervicksichtigen zu milssen, man nimmt nur die Sttitzwerte aus einer 

2 

bestimmten Umgebung von (A.. pJ, bei denen exp( - ad) grofier 
ist als eine vorgegebene Schranke. Die Werte der Koeffizienten wer- 
den dadurch nur unwesentlich verindert. Die LOsung des Gleichungs- 
systems (*=^0 sich durch die Bildung der Norma Igleichungen 


Q = 0 


rs 


in die Losung des linearen quadratischen Gleichungssystems 


T T 

A A c = A z 


umformen, wobei A eine n x 6 Matrix ist, c der Koeffizientenvektor 
und z der Vektor der Funktionswerte an den Sttitzstellen. Da die Sttitz- 
werte, die die ausgewerteten Interferogramme darstellen, zum Teil recht 

ungleichmSGig verteilt sind, ist es nicht ratsam, den Parameter a in 

2 

der Gewichtsfunktion W = exp(- a d ) fest vorzugeben, da sonst in einem 
Gebiet mit wenig Sttitzstellen u. U. eine zu starke Gewichtung der weiter 
entfernten Punkte, die ja jetzt zur Berechnung mit herangezogen werden 
mtissen, stattfindet. Der Parameter a wurde daher von der Sttitzpunkt- 
dichte abhtingig gemacht. 

Die Interferenzordnungsfunktion ist nun durch einen Satz von Polynom- 
koeffizienten an einem regelmafiigen Gitter gegeben und kann zur Er- 
zielung einer glatten Interpolationsfltiche leicht mit einem Splinever- 
fahren weiterverarbeitet werden [12]. 


Die Abb. 13ff zeigen als Anwendung der Polynomdarstellung im Teil b 
ein Linienbild eines Hologramms mit den eingezeichneten Gitterlinien 
und in Teil c die aus der Polynomdarstellung errechnete Interferenz- 
ordnungsfunktion in Hohenliniendarstellung. Vergleicht man beide Bilder, 
so stellt man fest, dafJ die Hohenlinien der Polynomdarstellung die 
Linien des Hologramms sehr gut reproduzieren und daB eine gewisse 
Giattung der Linien erfolgt ist. Durch die Wahl der Gittereinteilung 
fiir die Berechnung der Polynome ist ja eine kleinste auflosbare Struk- 
tur der FlcLche vorgegeben, wodurch eine TiefpaBfilterung der Inter - 
ferenzordnungsfunktion erfolgt, Diese Tiefpafifilterung ist durchaus 
erwiinscht und kann im Falle der Reifentesthologramme sogar fiir eine 
Reifenfehlererkennung ausgenutzt werden (s. Kap. 2.5), 

Eine weitere Eigenschaft dieser Darstellung ist in Abb. 13 u. 17 demonstriert. 
Hier konnte der Bereich des Hologramms, der durch einen im Sichtfeld 
liegenden Spreizer verdeckt war, nicht ausgewertet werden. Durch eine 
Anpassung der absoluten Numerierung des oberen und unteren Bildteils 
wurde durch die extrapolierende Eigenschaft der Polynomdarstellung 
die Interferenzordnungsfunktion auch in diesem verdeckten Bereich sehr 
gut beschrieben, was ein Vergleich mit dem Hologramm sofort bestatigt. 

Zur Berechnung der Polynomkoeffizienten wird das Programm PKF63 
benutzt, fhr die Darstellung von Schnittkurven der Polynomflache auf 
dem Sichtgerat das Programm PTST82 und zur Berechnung von geometri- 
schen Parametern, wie Kriimmungen und Gradienten der Polynomflache 
dient das Programm PTST84 (s. a. Programmbibliothek PKFLIB und 
SPLLIB in Teil II). 


Abb. 12: Beispielftir die Poly gonzugdarste Hung 

Interferenzstreifen 

a) am Interferogramm der Abb. 3a 

b) am Interferogramm der Abb 3b 
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2, 5 Summen und Differenzen von Interferogrammen 

Die in Kap. 2.4 beschriebene Polynomdarstellung der Interferenzordnungs- 
funktionen kann aufier fur die Berechnung von Zeitreihen bei den Stromungs- 
interferogrammen auch zu einer Mittelung iiber mehrere Interferogramme 
benutzt werden. 

Abb. 18 zeigt eine Anwendung auf Interferogramme eines turbulenten Gas- 
freistrahls. In Abb. 18b ist eine typische, mit dem Programm FAW 2 ein- 
gegebene Polygonzugdarstellung eines solchen Freistrahlinterferogramms 
dargestellt, wahrend die Abb. 18c eine iiber 20 Interferogramme gemittelte 
Interferenzordnungsfunktion zeigt. Man erkennt, da6 die turbulenten Schwank- 
kungen weitestgehend herausgemittelt sind und da(J die Funktion recht sym- 
metrisch erscheint, was auch theoretisch zu erwarten ist. 

In einer Anwendung der Polynomdarstellung auf die Fehlererkennung bei 
den Reifentesthologrammen wird die "TiefpaSeigenschaft" dieser Darstellung 
ausgenutzt, um kleinraumige Oberf lachendeformationen zu erkennen, durch 
die sich Reifenfehler hSufig darstellen. Zu diesem Zweck wind an den 
Punkten des Linienbildes die Differenz zu der berechneten PolynomflSche 
gebildet. Die Oberflachendeforraationen. die kleiner sind als die Maschen- 
weite der Polynomdarstellung werden sich durch grofiere Werte in diesem 
"Hochpafianteil" abzeichnen. Abb. 19 zeigt zwei solche Differenzbilder , 
wobei die Grofie der Abweichungen von der Polynomdarstellung durch unter- 
schiedliche Symbole kenntlich gemacht ist, 

Es zeigt sich, dafi sich an den "auff Slligen" SteUen auch tatsSchlich Reifen- 
fehler verbergen. Hier mud jedoch noch eine Strukturuntersuchung der 
Umgebung dieser Gebiete durchgefiihrt werden, um Aufschlud iiber die Art 
der Storung zu gewinnen und die Diagnose zu erhSrten (s.a. Kap. 2. 6). 



Abb. 18: 


Beispiel zur Interferogrammittelung 

a) Foto eines Interferogramms eines Uberschallfreistrahls 

b) mit dem Programm FAW2 von Hand eingegebene Strei- 
fenverlaufe des Interferogramms und Gitterlinien zur 
Berechnung der Polynomdarstellung 

c) liber 20 Interferogramme gemittelte Polynomdarstellung 
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Eine weitere Anwendung der Polynomdarstellung ist die Elimination 
von iiberlagerten Interferenzfeldern, Bei den Hologrammen ist durch 
die Art der Aufnahme und der Geometrie der Anordnung ein gewisses 
Interferenzmuster vorhanden, selbst wenn sich die ReifenoberflS.che 
nicht verformen wiirde. Eine Moglichkeit der Beriicksichtigung dieses 
Effektes wS.re, das Grundinterferenzmuster aus der Geometrie der 
Aufnahmeanordnung rechnerisch zu bestimmen und von den ausgewer- 
teten Interferenzordnungsfunktionen zu subtrahieren. Da jedoch nicht 
nur geometrische Parameter ftir die Form des Grundinterferenzmusters 
verantwortlich sind, sondern auch solche, die sich aus den Verformungs 
eigenschaften der Reifen ergeben und daher schwer zu Uberblicken sind, 
soil dieses sog. "Normhologramm" durch Mittelung iiber eine Reihe von 
Interferogrammen erhalten werden. Abb, 20 zeigt eine erste Approxi- 
mation eines solchen "Normhologramms" , das aus ftinf Interferogram- 
men gemittelt wurde. Die Abb. 20 zeigt aufierdem einige Differenzinter- 
ferogramme, die mit diesem Normhologramm gebildet wurden. 


2, 6 Erkennung spezieller Strukturen und Reifenfehler 

Nach Aussagen der Fa. Continental stellen sich im Hologramm die Reifen- 
fehler in folgender Weise dar: hohe Liniendichte, unregelmSBiger Linien- 
verlauf, ungleichma.(3ige Streifenbreiten, geschlossene Linien bzw. Linien- 
systeme relativ kleiner FlSche, sowie nicht auflosbare Strukturen, 

Zur Erkennung solcher Gebiete wurde das Programm REEL entwickelt, 
das aus der Polygonzugdarstellung des Hologramms bestimmte geometri- 
sche Parameter berechnet und mit Hilfe lokaler AuswahUcriterien die Stellen 
im Linienbild ermittelt, an denen Reifenfehler im obigen Sinne zu erwarten 
sind. Das Programm liefert zur Vorauswahl schwer auswertbarer Holo- 
gramme zunachst einige globale GroBen wie das VerhSltnis aus der Anzahl 
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geschlossener Linien zur Anzahl nicht geschlossener Linien, welches 
bei kontrastschwachen Hologrammen entsprechend gro(3e Werte annimmt . 

Die aufierdem berechnete mittlere Anzahl der Punkte, die zur Darstellung 
einer Langeneinheit der reduzierten Polygonziige benotigt wird, ist bei 
sehr zackigem oder gekrummtem Linienverlauf ebenfalls relativ grofi. Als 
weitere geometrische GroBen wurde fiir das gesamte Linienbild eine mittlere 
Liniendichte und ein mittlerer Linienabstand ermittelt, der zur Auswahl 
bestimmter Stellen im Bild dient, die fur Reifenfehler in Frage kommen. 

Eine Liste der geschlossenen Linien liefert einen Uberblick bzgl. Lage, 

Umfang, Flache und Form aller geschlossenen Polygonziige im Bild. Mit 

einem Parameter wird die Kreisahnlichkeit der Linien ausgedriickt, denii 

2 - 

der sogenannte Formfaktor berechnet sich gem. 4 n Fiache/(Umfang) 
und ist gleich 1 fiir die Kreisform. Mit Hilfe dieser geometrischen 
GroBe lassen sich auch die nach Anwendung des Programms BVSF (siehe 
Abschnitt 2.2.3) noch fehlerhaft verbundenen Linien ermitteln. 

Die fur die Reifenfehlererkennung wichtigen lokalen Storungen des Streifen- 
abstandes werden dadurch bestimmt, daB in jedem Bildpunkt normal zum 
jeweiligen mittleren Linienverlauf die Abstande zur linken und rechten 
Nachbarlinie berechnet und falls moglich der Quotient AQ dieser beiden 
GroBen gebildet werden (siehe Abb. 21a)). Berechnet man noch das Ver- 
haltnis AV der Abstandsquotienten zweier aufeinanderfolgender Punkte 
der Linien, so erhalt man bei entsprechend groBer Abweichung von 1 die 
Stellen mit starker Streifenbreitenvariation. Als Reifenfehler kommen 
aber nur diejenigen Orte in Frage, deren lokaler mittlerer Linienabstand 
(berechnet als Mittelwert der Abstande A1 , . . . , A4 in Abb. 21a.)) nicht 
groB gegenilber dem schon berechneten mittleren Linienabstand des Ge- 
samtbildes ist und die nicht in RandnShe Uegen. 
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Abb. 19: HochpaGanteil und Beispiele zur Strukturanalyse 

a) und c) Differenzen der PolynomflSche des Hologramms 
Nr. 16 bzw. Nr. 18 an den Linienpunkten des Streifenfeldes 
("HochpaBanteil") (die Abweichungen sind mit den Symbolen 

' . -+=0 6 markiert, wobei ' einer Differenz von 0. 075 

und # einer Differenz von 0. 75 Streifenbreiten entspricht) 

b) und d) nait Hi lie der Strukturuntersuchung erkannte Gebiete 
mit ileifenfehlern 










Abb. '20: Beispiel fiir Differenzinterferogramme 

a) Darstellung der Hohenlinien eines aus den Hologrammen 
der Abb. 13-17 gemittelten "Normhologramms" 

b) Hohenliniendarstellung der Differenz des Hologramms 
Nr. 16 zum Normhologramm a) 

c) wie b) Hologramm Nr. 18 d) wie b) Hologramm Nr. A4 



Gelegentlich stellen sich Reifenfehler als Zacken im Linienverlauf dar; 
deshalb haben wir diejenigen Stellen im Polygonzugbild ermittelt, an 
denen starke RichtungsSLnderungen der Linien zu verzeichnen sind. Die 
Zackenhohe ZH und Zackehbreite ZB wurden wie in Abb. 21b) darge- 
stellt, berechnet und mit Hilfe eines lokalen mittleren Linienabstandes 
(berechnet als Mittelwert der Abstande A1 , . . . , A4 in Abb. 21b)) nor- 
miert. Der Quotient ZH/ZB liefert eine Information liber die Form 
des Zackens, wShrend der Parameter ZP = RH ZH/ZB mit RH als 
relative ZackenhShe zur Auswahl bestimmter Zacken dient, die als 
Reifenfehler in Frage kommen. 



AQ.1 A1/A2 
AQ2 =A3/A4 

AV;:AQ2/AQ1 



Abb. 21 Skizzen zur Berechnung 

a) des Verhaitnisses AV und 

b) von ZackengroBen 



Es wurden in der Abb. 19 als Beispiele nur diejenigen Zacken als Kreise 
eingezeichnet, deren vom Programm REEL berechnete GroBen RH und 
ZP gewisse Schwellwerte iiberschreiten und die in einem Gebiet hoher 
Liniendichte liegen. Als Quadrate sind die Stellen mit starken Storungen 
im Streifenabstand gekennzeichnet, und wie man an den Ergebnissen sieht, 
werden die Reifenfehler , die sich in der Polygonzugdarstellung abzeichnen, 
recht gut erkannt. 


2. 7 Entzerrung von Objektgeometrien 

In vielen Fallen ist es zur Berechnung der Objektfunktion erforderlich, 
nach dem Vorliegen der Interferenzordnungsfunktion des Objektes eine 
Koordinatentransformation durchzufUhren. Bei den Stromungsinterfero- 
grammen gestaltet sich diese Transformation in der Regel einfacher, da 
hier nur eine Koordinate behandelt werden mufi., wobei die Interferenz- 
ordnung bei zweidimensionalen Stromungen durch eine lineare Beziehung 
und bei rotationssymmetrischen Stromungen durch die Losung einer 
Integralgleichung in die entsprechende Dichte umgerechnet wird. 

Bei den Reifentesthologrammen ist zusatzlich eine Entzerrung der Orts- 
koordinaten erforderlich, um die durch die Art der Aufnahme und Re- 
produktion der Hologramme sowie der Geometrie des Objektes erfolgte 
Transformation der Objektpunkte riickgangig zu machen und z.B. die 
Lokalisierung von bestimmten Strukturen in den Hologrammen in Reifen- 
koordinaten anzugeben. Abb. 22 zeigt eine Skizze der bei einem bestimmten 
Reifenformat auftretenden Verzerrungen, die durch die Zentralprojektion 
eines Torus in eine Ebene entstehen. Da eine Riicktransformation dieser 
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Abb, 22: Skizze zur Verzerrung der Objektkoordinaten durch die 

Ho logrammaufnahme 

^ ist der Winkel in Reifenlaufrichtung und a der Meri- 
dianwinkel in Grad 

( und a Sind hier fur die R eifendimension 12 R 22. 5 
berechnet) 



Abbildung nicht mathematisch geschlossen moglich ist, andererseits 
aber durch die Abweichungen des Reifens von der Torusform auch nur 
eine AnnS,herung sein wiirde, kbnnte man durch Ausmessen einiger auf 
demReifen markierter Punkte eine genauere Beriicksichtigung dieser 
Verzerrung erzielen. 

Die Durchfiihrung der Rticktransformation kann durch Interpolation auf 
den F IS-chen x (x.y) und y' (x,y), die durch die vorgegebenen Pafipunkte 
definiert sind, durch Anwendung von Splinealgorithmen [12] vorgenommen 
werden. Zur Berechnung einer Transformation der Ebene bei beliebig 
vorgegebenen PafSpunkten wurde das Programm KTRANF entwickelt 
(s. a, Teil II). 


2. 8 Rechenzeiten 

Die fiir die Auswertung eines Interferogramms benotigte Rechenzeit teilt 
sich auf in die Zeit zur Digitalisierung der Vorlagen, der Durchfiihrung 
eventuell notwendiger Vorverarbeitungsschritte, der Streifenextraktion 
und Numerierung der Linienbilder und schliefUich in die Zeit, die zur 
Weiterverarbeitung erforderlich ist, wie das Berechnen der Polynom- 
darstellung und die Fehleranalyse bei den Reifentesthologrammen. 

Die Digitalisiergeschwindigkeit hangt von der Abtastfrequenz, d. h. von 
der verwendeten Auflosung ab. Wegen des Zeilensprungverfahrens bei 
Videosystemen der CCIR-Norm mufJ nach der Digitalisierung einer Zeilen- 
folge 20 ms bzw. 40 ms gewartet werden, bis auf die nachsten Zeilen zu- 
gegriffen werden kann (s. Kap. 1. 3). Fiir ein Bild mit einer Auflosung von 
512 X 512 Punkten fallen daher etwa 10 Sekunden an, da jeweils zwei Zeilen 
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in einem Zuge digitalisiert werden. Wie bereits in Kap. 1. 3 angedeutet, 
konnten diese Zeiten auf ein Drittel reduziert werden, wenn in einem 
TV -Frame mehrere Zeilenfolgen abgetastet wiirden, was technisch mog- 
lich ware. Wenn auf Grund lichtschwacher Vorlagen zur Verbesserung des 
Signal-Rausch-Verhaltnisses eine zeUliche Signalmittelung durchgefiihrt 
werden mufi, was insbesondere bei den Hologrammen der Fall ist, multi- 
pliziert sich die Zeit mit der Anzahl der Mittelungen. Wiirde diese Signal- 
mittelung mit spezieller Hardware in Echtzeit ausgefiihrt (entsprechende 
Systeme sind in neuerer Zeit auf dem Markt), so wiirde lediglich eine 
Sekunde ftir eine Mittelung iiber 25 Bilder anfallen. 

Fiir die anderen in Kap. 2.1 besprochenen Bildvorverarbeitungsverfahren 
fallen bei einer BildgroBe von 512 x 512 Punkten ca. 30 Sek. Rechenzeit 
fiir die raumliche Mittelung (SM08), ca. 60 Sek. Rechenzeit fiir die Sub- 
traktion des mittleren Intensitatsprofils (SUBMF3) und ca. 60 Sek. fiir die 
Grauwerttransformation (^TREX 3) an. Diese Zeiten enthalten keine Zu- 
griffszeiten auf externe Speichermedien (I/O-Zeit), da bei einer Kombi- 
nation der Verfahren durch sequentieUe Abarbeitung eine Zwischenspeiche- 
rung der Daten nicht erforderlich ist. Lediglich die Anwendung des Pro- 
gramms SUBMF 3 erfordert ein Abspeichern der Daten. Hier wiirden dann 
noch ca. 20 Sek. I/O-Zeit hinzukommen. 

Falls zweidimensionale Filteroperationen (wie in Kap. 2.1.3 beschrieben) 
fiir die Bildvorverarbeitung angewendet werden, so fallen Rechenzeiten von 
ca. 2 0 Min. fiir eine Hin- und Riicktransformation mit Fast-Fourier- Trans-. 
form-Verfahren an, wenn die Berechnung auf dem PDP 11-Rechner durch- 
gefiihrt wird. Die Zeiten fiir die sich nun anschlieBenden Verarbeitungs- 
schritte wie Extraktion der Streifenrander , Numerierung und Polynom- 
darstellung sind von der Anzahl der Streifen im Interferenzbild abhangig 
und betragen fiir ein typisches Hologramm ca. 60 Sek. fiir die Streifen- 



extraktion und Redundanzreduktion der Polynomzuge (Kap. 2.2.2), ca. 

40 Sek. fiir die Numerierung (NUMl) und ca. 210 Sek. fiir die Berechnung 
der Po lynomflache (PKF63) (bei 117 Stiitzstellen wie z.B. in Abb. 13 ff). 

Die Reifenfehleranalyse mit dem Programm RFEL benotigt auch je nach 
Streifendichte unterschiedlich groOe Rechenzeiten, die in der GrdOenord- 
nung einiger Minuten liegen, wenn sie am PDPll -Rechner durchgefuhrt 
werden, weil hier sehr viele Operationen in REAL-Arithmetik implemen- 
tiert wurden.j um in der Entwicklungsphase des Programme keine zusatz- 
lichen Probleme durch Zeitoptimierung zu bekommen. 

Fiir die voUstandige Auswertung eines typischen R eifenfehlerhologramms 
mit zeitlicher Mittelung iiber 20 Bilder, raumlicher Mittelung, Streifen- 
extraktion, Berechnung der Polynomf lache , des Differenzhologramms 
ohne Reifenfehleranalyse fallen Rechenzeiten in der Grofienordnung 10-20 
Min. an. Fur die Auswertung eines Interferogramms ohne zeitliche und 
r&umliche Mittelung mit Berechnung der PolynomoberflS-che werden etwa 
5-10 Min. benotigt. Diese Zeiten sind fiir eine Serienauswertung recht 
hoch, weshalb wir uns fiir den Einsatz eines Arrayprocessors (AP) am 
System entschieden haben. Insbesondere bei der Bildvorverarbeitung und' 
der Berechnung der Polynomoberf lache konnen die Zeiten nach Einsatz 
des Arrayprocessors um den Faktor 50-100 reduziert werden, wie Vergleichs- 
rechnungen an einem AP gezeigt haben. Auch die zeitintensiven Segmente 
der anderen Programme sollen auf den AP verlagert werden, wodurch 
wir eine Reduzierung der Gesamtauswertezeit auf ca. 20-40 Sek. pro 
Interferogramm bzw. Hologramm erwarten. Durch den Einsatz des AP 
sind auch Verfahren der zweidimensionalen Filterung zur Bildvorverarbei- 
tung in der Serienauswertung realisierbar . 



3. 


Programme zur Interferogrammauswertung 


Die im vorausgegangenen Text bei der Besprechung der einzelnen Ver- 
fahrensschritte zitierten Programme sind in einem zweiten Teil dieses 
Berichtes zusammengefafit [27]. In die Programmlisten wurde zum 
grofien Teil erlSuternder Kommentar eingefugt, der ein besseres Ver- 
stSndnis der Programme ermoglicht. Dem Band ist aufierdem ein alpha- 
betisches Verzeichnis der Programme nachgestellt, das ein bequemes 
Aufsuchen bestimmter Programme ermoglicht. 

Urspriinglich war im Rahmen des Forschungsprojektes vorgesehen, das 
Programmsystem in ein entsprechendes Auswertesystem beim Projekt- 
partner Continental Gummiwerke zu integrieren. Weil jedoch die. ent- 
sprechenden Investitionen dort nicht durchgefiihrt wurden, konnte dieser 
Teil des Projektes nicht realisiert werden. 
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